/*BM字符串匹配算法*/
/*code by CG cg45.cn
* ACM yctc
*2008 12 20
*/
#include <stdio.h>
#include "string.h"
#include "stdlib.h"
 
#define LEN 256
/*LEN 使用一个character set ASCII编码使用一个字节表示字符，UNICODE要另
*外考虑，BM算法优势所在*/
 
/*BM() BM算法基本功能函数
*输入:
*	char *s 匹配串
*	char *p 模式串
*	int index 匹配开始索引
*	int post[] 辅助数组
*返回:
*	int 下一个匹配开始的索引，匹配失败返回-1
*/
int BM(char *s, char *p, int index, int post[]) {
   int len = strlen(s);
   int i,j, next;
   i = strlen(p)-1;/*字符串长度减1*/
   j = index+strlen(p)-1;
/*第一次调用 BM() 时 index = 0，因
*为下面的 for 循环是从模式串的末尾开始比较，所以匹配串的初始比较位
*置应该是从开头数模式串长度个位置开始。*/
   for(; i>=0; i--, j--) {
      if(s[j] != p[i]){/*第一个字符的匹配*/
        break;
    }
   }/*for*/
 
   if(i<0) /*匹配完毕?*/
     return 0; /*匹配成功*/
 
   else if(post[s[j]]>0 && post[s[j]] < i) {
/*当出现不匹配时，查看匹配串当前位置的字符有没有出现在模式串中*/
     next = index + i - post[s[j]];
 
/*index 是当前的匹配串起始偏移量，i 是模式串还剩的比较字串数目，
* post[s[j]]是所出现的第一个不匹配的字符在匹配串中的位置。
*这样下次比较就从匹配串中出现 s[j] 的位置开始比较
*/
   }
   else next = index + 1;
 
   if(next > len-strlen(p))
     return -1; /*匹配失败，无法进行下一次匹配*/
   else
     return next; /*匹配失败，需要下一次匹配*/
 }/*BM*/
 
 /*测试,匹配串 和 模式串都使用小写字符*/
 int main()
 {
 int post[LEN]={0}; /*辅助数组=字符集大小*/

 char *src="aaaabbbaababaaabbabb";/*测试字符串*/
 char *patten="aabbabb";

 int i, next, index, pos=0;/*初始化索引标志*/

 for(i=0; i<strlen(patten); i++) /*构造辅助数组*/
    post[patten[i]]=i;

 index = BM(src, patten, 0, post);/*第一次匹配，从0位置开始,获得NEXT*/

 while(!(index == -1 || index == 0)) /*循环直到匹配成功*/
 {
   next = index;
   index = BM(src, patten, next, post);/*下一次BM匹配*/
 }

 if(index == -1){ /*faild*/
    printf("Match faild\n");
 }

 if(index == 0){ /* OK */
    printf("the index is: %d.\n", next);
 }
 return 0;
 }
